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8 : Abstract 

| In this article we study the general structure and special properties of the Schwinger-Dyson equa- 

■ tion for the gluon propagator constructed with the pinch technique, together with the question of 

how to obtain infrared finite solutions, associated with the generation of an effective gluon mass. 

Exploiting the known all-order correspondence between the pinch technique and the background 

<^ ' field method, we demonstrate that, contrary to the standard formulation, the non-perturbative 

O ' 

gluon self-energy is transverse order-by-order in the dressed loop expansion, and separately for 
\Q ■ 

gluonic and ghost contributions. We next present a comprehensive review of several subtle is- 



Oh 

CD 



Qh| sues relevant to the search of infrared finite solutions, paying particular attention to the role of 

the seagull graph in enforcing transversality, the necessity of introducing massless poles in the 
three-gluon vertex, and the incorporation of the correct renormalization group properties. In ad- 
^ ; dition, we present a method for regulating the seagull-type contributions based on dimensional 

regularization; its applicability depends crucially on the asymptotic behavior of the solutions in 
the deep ultraviolet, and in particular on the anomalous dimension of the dynamically generated 
gluon mass. A linearized version of the truncated Schwinger-Dyson equation is derived, using a 
vertex that satisfies the required Ward identity and contains massless poles belonging to different 
Lorentz structures. The resulting integral equation is then solved numerically, the infrared and 
ultraviolet properties of the obtained solutions are examined in detail, and the allowed range for the 
effective gluon mass is determined. Various open questions and possible connections with different 
approaches in the literature are discussed. 

PACS numbers: 12.38.Lg, 12.38.Aw 
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I. INTRODUCTION 



The generation of mass gaps in QCD is one of the most fundamental problems in particle 
physics. In part the difficulty lies in the fact that the symmetries governing the QCD 
Lagrangian prohibit the appearance of mass terms for all fundamental degrees of freedom 
at tree-level and, provided that these symmetries are not violated through the procedure of 
regularization, this masslessness persists to all orders in perturbation theory. Thus, mass 
generation in QCD becomes an inherently non-perturbative problem, whose tackling requires 
the employment of rather sophisticated calculational tools and approximation schemes Q]. 

Whereas the generation of quark masses is intimately connected with the breaking of 
chiral symmetry |2j, it was argued long ago that the non-perturbative QCD dynamics lead 
to the generation of an effective gluon mass, while the local gauge invariance of the theory 
remains intact . This gluon "mass' is not a directly measurable quantity, 

but must be related to other physical parameters such as the string tension, glueball masses, 
or the QCD vacuum energy lflj . and furnishes, at least in principle, a regulator for all 
infrared (IR) divergences of QCD. 

The concept of a dynamically generated gluon mass, its field theoretic realization, and a 
plethora of physical and technical issues associated with it, have been explored in great detail 
in a classic paper was the insistence 

on preserving, at every level of approximation, crucial properties such as gauge-invariance, 
gauge-independence, and invariance under the renormalization-group (RG). With this mo- 
tivation, an effective gluon propagator, A^, was derived through the systematic rearrange- 
ment of Feynman graphs, a procedure that is now known in the literature as the "pinch 



technique" (PT) (l2j, 3, Il4j . The self-energy IT^, of this propagator is gauge- independent 



and captures the leading logarithms of the theory, exactly as happens with the vacuum po- 
larization in QED. The central result of was that, when solving the Schwinger- Dyson 
(SD) equation governing the PT propagator, and under special assumptions for the form 
of the three-gluon vertex, one finds solutions that are free of the Landau singularity, and 
reach a finite value in the deep IR. These solutions may be successfully fitted by a massive 
propagator, with the crucial characteristic, encoded in the corresponding SD equation, that 
the mass employed is not "hard", but depends non-trivially on the momentum transfer, 
vanishing sufficiently fast in the deep ultraviolet (UV). From the dimensionfull massive so- 
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lutions one may define a dimensionless quantity, which constitutes the generalization in a 
non-Abelian context of the universal (process-independent) QED effective charge. The QCD 
effective charge so obtained displays asymptotic freedom in the UV, whereas in the IR it 
"freezes" at a finite value. 

Various independent field theoretic studies 
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18, 



19 



20 



21(, spanning over a 



quarter of a century a 



computations 



22, 2c 



so corroborate some type of gluon mass generation. In addition, lattice 
24| reveal the onset of non-perturbative effects, which in principle can 
be modelled by means of effectively massive gluon propagators. It is important to emphasize 
that the massive gluon propagator derived in |ll[ describes successfully nucleon-nucleon 
scattering when inserted, rather heuristically, into the two-gluon exchange model j^; for 
additional phenomenological applications, see l26l. Furthermore, several theoretical studies 



based on a-priori very distinct approaches 



29 



30, 



31 



32 



33, 



34, 



35, 



37 



support the notion of the "freezing" of the QCD running coupling in the deep IR (but do 
not agree, in general, on its actual value). 

In recent years there has been significant progress in our understanding of the PT con- 
struction in general jsjj, and the properties of the resulting effective Green's functions in 
particular. The extension of the PT to all orders was carried out in and the known 



4J,|43J was shown to persist to all 



one- [40J and two-loop [41] connection with the BFM 
orders. From the practical point of view the established connection permits the direct cal- 
culation with a set of concrete Feynman rules, and enables one to prove all-order results, 
exploiting the powerful formal machinery of the BFM. In what follows we will refer to the 
framework emerging from the synergy between PT and BFM as the "PT-BFM scheme" . 

The aim of this article is threefold: First, we initiate a systematic treatment of the SD 
equations within the PT-BFM scheme, with particular emphasis on the manifestly gauge- 
invariant truncation it offers. Second, we discuss various field-theoretic issues relevant to 
the study of gluon mass generation in the context of SD equations in general. Third, we 
analyze in detail the SD equation obtained as the first non-trivial approximation in the 
aforementioned truncation scheme, and search for infrared finite solutions. 

Regarding our first objective, let us point out that one of the most distinct features of 
the PT-BFM scheme is the special way in which the transversality of the background gluon 
self-energy is realized. In particular, the study of the non-perturbative, SD-type of 
equation obeyed by reveals that, by virtue of the Abelian-like Ward Identities (WI) 
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satisfied by the vertices involved, the transversality is preserved without the inclusion of 
ghosts. Put in another way, gluonic and ghost contributions are separately transverse. In 
addition, transversality is enforced without mixing the orders in the usual "dressed-loop" ex- 
pansion: the "one-loop-dressed" and "two-loop-dressed" sets of diagrams are independently 
transverse. This is to be contrasted to what happens in the usual gauge-fixing scheme 
of the covariant renormalizable gauges, where the inclusion of the ghost is crucial for the 
transversality already at the level of the one-loop perturbative calculation. This particular 
transversality property of the BFM self-energy is known at the level of the one-loop calcu- 
lation however, to the best of our knowledge, its all-order generalization presented in 
Sec. in] appears for the first time in the literature. The importance of this property in the 
context of SD equation is that it allows for a meaningful first approximation: instead of 
the system of coupled equations involving gluon and ghost propagators, one may consider 
only the subset containing gluons, without compromising the crucial property of transver- 
sality. More generally, one can envisage a systematic dressed loop expansion, maintaining 
transversality manifest at every level of approximation. 

Instrumental for some of the developments mentioned above has been a set of non-trivial 
identities j^J, relating the BFM n-point functions to the corresponding conventional n-point 
functions in the covariant renormalizable gauges, to all orders in perturbation theory. These 
identities, to be referred to as Background-Quantum identities (BQIs) |39j |. are expected 
to play a fundamental role in addressing one of the most central issues in the context of 
the PT-BFM scheme, namely the actual construction of a new SD series. Specifically, as 



is known already from the two-loop analysis (41 1, the PT-BFM gluon self-energy 14^ is 
expressed in terms of Feynman diagrams containing the conventional gluon self-energy 14^; 
this fact is generic, as the all-order diagrammatic representation of IT^ demonstrates (see 
Sec. |nj). Clearly, in order to arrive at a genuine SD equation for IT^, one must carry out 
the substitution — > Yi^ v inside the loops. It is still an open question whether such a 
replacement can be implemented self-consistently to all orders; a preliminary view of how 
this might work out is presented in Sec. ITT1 

Turning to the analysis of the SD equations and the search for infrared finite solutions, 
after setting up the appropriate theoretical stage in Sec. IIII1 in the next two sections we 
eventually study a linearized version of the equation governing the PT propagator V, in 
the spirit of |TJj. Although several of the techniques developed there are adopted virtu- 
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ally unchanged in the present work, there are important theoretical and phenomenological 
differences, which we summarize below. 

(i) The role of the ghosts: As has become clear from the detailed study of the correspon- 
dence between PT and BFM, the rearrangement of graphs (or sets of graphs ) im plemented 
by the PT generates dynamically the characteristic ghost sector of the BFM 14 1 . However, 
Ice the o ^al oee-.oe, derilon of t he PT S e lt -ene rgy Q th e c*J». were car- 
ried out in the context of the ghost-free light cone gauge, the distinction between gluonic 
and ghost contributions was not so obvious. As a result, in the heuristic derivation of the 
corresponding SD equation, all contributions were treated as gluonic. This is reflected in the 
fact that the coefficient multiplying the characteristic term ln(g 2 +4m 2 (g 2 )) appearing in the 
solutions (e.g. the standard RG logarithm supplemented with the non-perturbative mass) is 
precisely b = HCa/^Stt 2 , namely the coefficient of the one-loop QCD (3 function. Instead, 
the PT-BFM correspondence reveals that the purely gluonic contributions is b = 10Ca/487t 2 . 
Needless to say, the point is not so much the minor numerical discrepancy in the coefficients 
multiplying the logs, but rather the possibility that the ghost dynamics may behave in a 
completely different way in the IR. Thus, whilst the ghosts will eventually furnish the miss- 
ing Cyi/487r 2 asymptotically, their IR contribution may deviate from the massive logarithm 
given above, inducing qualitative changes in the form of the full gluon self-energy. Reversing 
the argument, in order to actually obtain solutions of the type 61n(g 2 + 4m 2 (g 2 )) from the 
coupled gluon-ghost system of SD equations, a very delicate interplay between gluons and 
ghosts must take place. 

(ii) Form of the three-gluon vertex: It is well-known that in order to obtain dynamically 
generated masses one must allow for the presence of massless poles in the corresponding ex- 



11] 



pression for the three-gluon vertex |45[ • The effective vertex used in the SD equation of 
was the bare three-gluon vertex obtained from the Lagrangian of the (non-renormalizable) 

nn 

massive gauge- invariant Yang-Mills model |4fiL 14 7| ; it contains kinematic poles, whose di- 
mensionality is partially compensated by the explicit appearance of a hard mass term in 
the numerator. Instead, we use a gauge-technique inspired Ansatz for the vertex [3] , which 
also contains kinematic poles, but their dimensionality is saturated solely by appropriate 
combinations of the momenta involved, with no explicit reference to mass terms, thus being 
closer to what one might expect to obtain within QCD. We hasten to emphasize that our 
Ansatz for the vertex is completely phenomenological, and is not derived from any dynam- 
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ical principle, other than the WI that it satisfies, nor does it exhaust the possible Lorentz 
structures. What we hope to obtain by resorting to such a simplified vertex is a manageable 
SD equation, that will allow us to study in detail the complicated interplay of the various 
components, and get a feel for the dependence of the solutions on the form of the vertex 
used. 

(iii) Seagull regularization: As was explained in [11] , the integral equation describing gluon 
mass generation is supplemented by a non-trivial constraint, expressing A _1 (0) in terms of 
(quadratically divergent) seagull-like contributions; after its regularization, this constraint 
will restrict severely the number of possible solutions. The phenomenological "glueball 
regularization" employed in |h| was based on an elaborate connection between the seagull 
contributions, the massive Yang-Mills model, and the finite vacuum expectation value of 
a scalar field creating glueball states. Instead, the regulation we introduce in this work 
is based solely on dimensional regularization. In particular, the non-perturbative seagull 
contributions are regulated by subtracting from them the elementary integral J d 4 k/k 2 , 
which vanishes in dimensional regularization. It turns out that this subtraction is sufficient 
to regulate the expression for A _1 (0), provided that the momentum-dependent mass vanishes 
"sufficiently fast" in the deep UV. In turn, this required asymptotic behavior restricts the 
values of the parameters appearing in the integral equation. 

(iv) Type of solutions: For relatively moderate values of A _1 (0), which at the level of 
the integral equation is treated as an input, the type of solutions emerging may be fitted 
with great accuracy by means of a monotonically decreasing dynamical mass and a running 
coupling, exactly as advocated in jllj] . However, as one decreases A _1 (0) beyond a critical 
value, a new class of qualitatively different solutions begins to emerge. These solutions 
are also finite in the entire range of momentum, but they display a sharp increase in the 
deep IR, and the corresponding plateau-like range, associated with the "freezing" , becomes 
increasingly narrower. 

The paper is organized as follows: In Sec. |H] we review the BT-BFM scheme, and the 
structure of the all-order gluon self-energy. In Sec. IIHI we present some general consider- 
ations pertinent to the search for IR finite solutions, with particular emphasis on the role 
of transversality, the kinematic poles in the vertex employed, the restoration of the correct 
UV behavior at the level of the SD equation, and the regularization of the seagull terms. In 
Sec. II VI we derive the linearized integral equation and discuss in detail several of its charac- 
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teristics. Sec. IS contains the numerical analysis, focusing particularly on the appearance of 
two types of solutions, as mentioned above. Finally, in Sec. I VII we discuss connection of this 
work with other approaches in the literature, outline various possible future directions, and 
summarize our conclusions. 

II. THE PT-BFM SCHEME 

In this section we study the structure of the effective gluon self-energy obtained within 
the PT-BFM framework. In the first subsection we present a brief overview of the PT and 
its connection with the BFM. The discussion presented here is meant to serve as a brief 
reminder; for a more complete treatment the reader is referred to the extensive literature 
on the subject. In the second subsection we present the all-order diagrammatic structure 
of the gluon propagator. In the third subsection we first derive an elementary WI, valid in 
the ghost sector of the BFM, and then demonstrate that, to all orders, the contributions of 
gluonic and ghost loops to the effective gluon self-energy are separately transverse. Finally, 
in the last subsection we present a preliminary view of how the PT may eventually lead to 
a new SD series. 



A. The connection between PT and BFM 



The PT 
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12j is a well-defined algorithm that exploits systematically the symmetries 
built into physical observables, such as S'-matrix elements, in order to construct new, ef- 
fective Green's functions endowed with very special properties. Most importantly, they are 
independent of the gauge-fixing parameter, and satisfy naive (ghost-free, QED-like) Wis 
instead of the usual Slavnov- Taylor identities. The basic observation, which essentially de- 
fines the PT, is that there exists a fundamental cancellation between sets of diagrams with 
different kinematic properties, such as self-energies, vertices, and boxes. This cancellation is 

n 

driven by the underlying BRST symmetry |42| , and is triggered when a very particular sub- 
set of the longitudinal momenta circulating inside vertex and box diagrams generate out of 
them (by "pinching" out internal lines) propagator-like terms. The latter are reassigned to 
conventional self-energy graphs, in order to give rise to the aforementioned effective Green's 
functions. 
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The longitudinal momenta responsible for these diagrammatic rearrangements stem either 
(a) from the bare gluon propagators contained inside the various Feynman diagrams, 

kitkij 



v< - (i - 



(2.1) 



and/or (b) from the "pinching part" (g,pi,p 2 ) appearin g in the characteristic decom 



position of the bare three-gluon vertex T 



eab 



9f eab ^au V into 



T af iu(q,Pi,P2) = r^(g,pi,p 2 ) +Tl fMU (q,p 1 ,p 2 ). 

r ^U(9,Pl,P2) = (Pi -P2)a9p, + 2qu9av ~ ^9o 
F L^(<?,Pl,P2) = Plvga^-Px^oa,- 



(2.2) 



The case of the gluon self-energy is of particular interest. Defining the transverse projector 



q 

we have for the full gluon propagator in the Feynman gauge 



P^(g)A(g 2 ) + 



The scalar function A(q 2 ) is related to the all-order gluon self-energy U^ u (q), 



(2-3) 



(2.4) 



IW<?) = P^(g)n(g 2 



(2.5) 



through 



A(g 2 



(2.6) 



q 2 + iU(q 2 ) 

Notice that the way H IJ , 1/ (q) has been defined in ()2.fij) (e.g. with the imaginary factor i in 
front), it is given simply by the corresponding Feynman diagrams in Minkowski space. The 
inverse of the full gluon propagator has the form 



(2.7) 



or, equivalently, 

A'M^ig^-U^q). (2.8) 

The PT construction of the effective one-loop self-energy Il\ll(q) can be most easily con- 
structed directly in the Feynman gauge. It amounts to adding to the conventional one-loop 
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n$(9)= \ 




+ 



+ 



WW! 



WW* 



(b) 



!ww 



+ 2 




(c) 



FIG. 1: The PT self-energy at one-loop. 

n[!J/(g) ( Fig. H (a) and (b)) the pinch contributions coming from vertex graphs, shown 
schematically in (c). Then, the final result is 

,2 ' 



Q 



(2.9) 



where b = 11Ca/487t 2 is the first coefficient of the QCD /3-function. 

Evidently, due to the Abelian Wis satisfied by the PT effective Green's functions, the 
new propagator-like quantity A _1 (g 2 ) absorbs all the RG-logs, exactly as happens in QED 
with the photon self-energy. Equivalently, since Z g and Za, the renormalization constants of 
the gauge-coupling and the effective self-energy, respectively, satisfy the QED relation Z„ = 
Z^ 1//2 , the product d(q 2 ) = g 2 A(q 2 ) forms a RG-invariant (//-independent) quantity |50l l51|: 
for large momenta q 2 , 



d(q 2 



q 2 



where g 2 {q 2 ) is the RG-invariant effective charge of QCD, 

^ 2 1 



3=V 



9 



(2.10) 



(2.11) 



1 + bg 2 In (q 2 /n 2 ) b In (q 2 /A 2 ) ' 

Of central importance for what follows is the connection between the PT and the BFM. 
The latter is a special gauge-fixing procedure, implemented at the level of the generat- 
ing functional. In particular, it preserves the symmetry of the action under ordinary 
gauge transformations with respect to the background (classical) gauge field A^, while 
the quantum gauge fields A ^ appearing in the loops transform homogeneously under the 
gauge group, i.e., as ordinary matter fields which happened to be assigned to the adjoint 
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representation 



As a result of the background gauge symmetry, the n-point func- 



tions (0\T 



(xi)A a ^ 2 (X2) ■ ■ ■ A a ™ n (x n ) |0) are gauge- invariant, in the sense that they satisfy 
naive, QED-like Wis. Notice, however, that they are not gauge-independent, because they 
depend explicitly on the quantum gauge-fixing parameter £q used to define the tree-level 



propagators of the quantum gluons and 
and two background gluons, respectively 



;he three- and four-gluon vertices involving one 
43I ]. The connection between PT and BFM may 



be stated as follows: The (gauge-independent) PT effective n-point functions (n = 2,3,4) 
coincide with the (gauge-dependent) BFM n-point functions (n = 2,3,4 background gluons 
entering) provided that the latter are computed at £q = 1 (e.g. setting £q = 1 in the 
Feynman rules of the Appendix). This connection was first established at one- loop level j^j], 
and was recently shown to persist to all orders in perturbation theory [l^ . 



B. The SD equation of the effective gluon self-energy 

The structure of the effective gluon self-energy, as it emerges from the all-order PT- 
BFM correspondence, can be written in a closed non-perturbative form, which coincides 
with the SD equation for A, derived formally from the BFM path integral using functional 



techniques 53]. 



In what follows we assume dimensional regularization, and employ the short-hand nota- 
tion [dk] = d d k/ (2n) d , where d = 4 — e is the dimension of space-time. We refer to diagrams 
containing one explicit integration over virtual momenta as "one-loop dressed" and those 
with two integrations as "two-loop dressed". 

We will classify the corresponding diagrams into four categories: one-loop dressed gluonic 
contribution (group a), one- loop dressed ghost contribution (group b), two-loop dressed 
gluonic contribution (group c), and two- loop dressed ghost contribution (group d). 

The closed expressions corresponding to the two diagrams of Fig. El are given by 



1 f ~ ~hf'r' 
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FIG. 2: The gluonic contribution at one-loop dressed expansion. Wavy lines with grey blobs 
represent full-quantum gluon propagators. All external wavy lines (ending with a vertical line) 
are background gluons. The black dots are the tree-level vertices in the BFM, and the white-blob 
denote the full three-gluon vertex with one background gluon. 




FIG. 3: The ghost sector at one-loop dressed expansion. Dashed lines with grey blobs denote 
full- ghost propagators, while the white blob represents the full background gluon- ghost vertex 

For the Fig. |3] we have 

^ f ~be'x l 

n ?^)| bl = - J[dk)r^D ee ,(k)T v D xx ,(k + q), 
fi»| b2 = - f[dk\Y^D ex {k). (2.13) 
The two-loop dressed gluonic contribution, Fig. HI reads 

n£(?)| ci - -/ J[dk][d£]r^A%\k^ 

K b M\ C2 = \ J WAS'(^ + ^^A:;(i)f ^A^(f + g ) . 

(2.14) 

11 



(ca) 



FIG. 4: Purely gluonic graphs relevant for the two-loop dressed expansion. The black blob rep- 
resents the full conventional three-gluon vertex, while the white blobs denote three or four-gluon 
vertices with one external background leg. 



l + k 




FIG. 5: The ghost sector contribution to the two-loop dressed expansion. The black blobs are 
the conventional full gluon-ghost vertices, while the white ones represent ghost vertices with an 
(external) background gluon and two ghosts. 

The last group represents the two-loop dressed ghost contribution, Fig. EJ and is written 
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as 

K (?)| dl = - j JWM T 7a eXD cAk + f)A» (k)D xx ,(£ - g)r ra , , 

nj(<?)| d2 = - / y^«f^^ 
nj(?)| d3 = - / ^A«f^ 

nj(g)| d4 = - | ^p^^fc + ^A^ (2.15) 

Notice that, (i) as explained in the Introduction, the propagators appearing inside the 
loops are quantum ones, and (ii) there are two general types of vertices, those where all 
incoming fields are quantum, and those where one of the incoming fields is background. 



C. Special transversality properties 

It is well-known that in the conventional formulation, the diagram containing the ghost- 
loop (graph (b) in Fig.QJ is instrumental for the transversality of (q). On the other hand, 
in the PT-BFM scheme, due to the special Feynman rules (see Appendix), the contributions 
of graphs (a) and (b) are individually transverse. Specifically, keeping only the logarithmic 
terms, one has |___| 

n<5(?) = m g) p„(«) . n*»( 9 ) = ^ m g) p„(,) . (2.16) 

In this subsection we will show that, by virtue of the all-order WI satisfied by the full 
vertices appearing in the diagrams defining H lll/ (q), Figs.(J____J), the above property is valid 
non-perturbatively, and that gluonic and ghost contributions are separately transverse. In 
addition, the one-loop and two-loop dressed diagrams do not mix. This is to be contrasted 
to what happens to the conventional case, where the orders of the loop expansion also mix. 

There are four fully dressed vertices with one incoming background gluon appearing in the 

_____ — abc — acb — abed ~ cdba 

diagrammatic definition of !!„,,,( g), in Figs.(J_0|5j): IT T , IT «, IT . As is known from 



formal considerations (see 



and references therein), the WI obtained when contracting 



such vertices with the momentum carried by the background gluon retain to all-orders the 

~ abc — acb 

same form as at tree-level. The tree-level WI for IT 1 ^ and are simply 

q^%(q, Pl ,p 2 ) = 9 r bc (pl - pl)g aP , q^f(q, Pl ,p 2 ) = gf*°{& - p\) . (2.17) 
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In addition, the fact that the four-gluon vertex with one incoming background gluon, and 
the conventional one (four quantum gluons) coincide at tree-level, as shown in Eq. ([A-5j) . 



furnishes the corresponding tree-level WI for F n (see, for instance [54|). Therefore, the 



fii/af3 

— cdba 

only ingredient missing is the tree-level identity satisfied by F ; we now proceed to its 
derivation. 

Contracting the bare vertex T cdba = —ig 2 f acx f xdb g fl v, shown in Fig. IA-91 with the mo- 
mentum ft carried by the background gluon, we have that 

tff ?°(<zi, ft, g 3 , ft) = -ig 2 f acx f xdb qiv = ig 2 (f abx f cdx + f°**f»>*) qiv , (2.18) 

where we have used the Jaccobi identity 

jabx jcdx _|_ jacx jdbx _|_ jadx jbcx q ^2 

Next we use that 

= _ ^abxjcdx + jacxjdhx + jadxjbcx^ ^ + g ^ ^ ^ ^ 

and add it by parts to (|2.18|) . obtaining 

„fij^cdba -2 \ rcdx faxb„ , ecbx iiadx , rcax fxdbt , \ ] 

Ql L tiu — W [J J QAu + J J qAu + J J Ift + ftjz/j 

= -ig [f cdx T a u xb (q A , q 2 + ft, q 3 ) + f cbx T adx {q^ q 2 , q 3 + ft) + f cax T x u db (q A + ft, q 2 , q 3 )] . 

(2.21) 

Armed with the above tree-level results, we proceed to state the four fundamental all- 
order Wis. First, the WI of the three-field vertices, where on the RHS we have differences 
of inverse propagators, are given by 

tfO^**) = gf abc [A-J(q 2 ) - A-J(q 3 )} , 

q^\q 2 , ft, q 3 ) = gf abc [D~\q 2 ) - D" 1 ^)] . (2.22) 

Then, the WI of the four-field vertices, where on the RHS we have sums of three trilinear 
vertices, with appropriately shifted arguments, are 

?i^i(9i> 92, ft, ft) = wf abx K%(^ 94, ft + ft) 

+ tgf acx T db u x Q (q 4 ,q 2j q 1 + q 3 ) 



+ igf adx T b ^{q 2 ,q 3 ,q 1 + q i ), (2.23) 
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and 

gflOgi, q 2 , 93, 94) = -igf cdx T a u xb (q 4 , q 2 + qi , q 3 ) 

^gr bx T a u d *(q A ,q 2 ,q 3 + qi ) 
-igf cax r^(q 4 + qi ,q 2 ,q 3 ). (2.24) 

Notice that Eq. (j2.24|) is the all-order generalization of (j2.21|) . 

With the above WI we can prove that the four groups presented before are independently 
transverse. We start with group (a): we contract graph (ai) using the first all-order WI 
of (J2.22j) . whereas for (a 2 ) we simply compute the divergence of the tree- level vertex (jA-6|) . 
arriving at 

g^t(g)| ai = C A g 2 5 ab q, J[dk] A' p (k) , 

Q u nf»(Q)L 2 = -C A g 2 5 ab q, J[dk) A» p (k) , (2.25) 

Thus, 

^(nJ(g)| ai + n^(g)| a2 )=o. (2.26) 

Similarly, the one-loop-dressed ghost contributions of group (b) give upon contraction 

9"n£(?)| bl = 2C A g 2 5 ab q u J[dk]D{k), 

q v UfM\ b2 = -2C A g 2 6 ab q u J[dk}D(k). (2.27) 

and so 

^(^)| bl +n^(g)| b2 )=0. (2.28) 

The two-loop dressed demonstration is slightly more involved, but essentially straightfor- 
ward. We begin with the two-loop gluonic contributions (group c). The action of q v on the 

-bx'e'c' 

all-order four-gluon vertex F , (3 , Q ,,(— q, £ + q, — k — £, k) appearing on the RHS of the first 
equation in f)2. 14j) may be obtained from ()2.23jl . through the following definition of momenta, 
qi = —q, q 2 = £ + q, q 3 = —k — £, g 4 = k, and corresponding relabellings of Lorentz and 
color indices. Then, 

q v fifM\ cl = \w j j[dk][d£}f^A aa \k)A^\k + £)A^(£ + q)[ 

f b ™T e ^,(-k - £, k, £) + f bei T% l0 ,(k, £ + q ,-k-£-q) + 
f b ™T^, a ,(£ + q,-k-£,k-q)]. 

(2.29) 
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It is not difficult to verify that, after judicious shifting and relabelling of the integration 
momenta and of the "dummy" Lorentz and color indices, the three terms on the RHS of 
(|2.29|) are in fact equal. Thus, 

q"U%(q)\ ci = l -igf™ j ' f[dk][dl]Y^A aa \k)A^ ' (k + £)A^' (£ + q)Tf alj ,(-k - £ : k,£) . 

(2.30) 

The other graph gives 

q v UfM\ C2 = \^f bxt jj[dk}[di]f^ A^\k)A^\k + £)T% a ,(£,-k-£,k)A^(£ + q) 

- \wf xi J J[dk}[d£]f^ A^{k)A^\k + £)T% a ,{£-k-£,k)A^{£). 

(2.31) 

Evidently the second term on the RHS of ()2.31|) vanishes identically, since the integral is 
independent of q, and therefore the free Lorentz index fi cannot be saturated. 

We then observe that, due to the full Bose symmetry of the conventional three-gluon 
vertex, Tp? a , ,(-k - £, k, £) = Fy^, Q ,(£, -k -£,k), and therefore, finally, 

9"( fi St(5)U + fi J(9)lo a )=°- ( 2 ' 32 ) 

Finally, we turn to the two-loop-dressed ghost-graphs (group d). For the calculation of 
the divergence of graph (d x ) we use Eq. fl2.24| ): this WI generates three distinct terms, i.e. 

<?Xt(g)| di = -i9 ff[dk][de]f a ^D(k + e)A^(k)D(£-q)[ 

f eM T xic( q _^_ k _ q ^ k + i) + f**FF(q -£,-k,k + e-q) 

+ fbxi T iec(_^ k + ^ _ 

(2.33) 

Each one of these three terms can be easily shown to cancel exactly against the individual 
divergences of the remaining three graphs. To see this in detail, use the first WI of ()2.22|) 
for graph (d 2 ), and the second of (|2.22|) for graphs (d 3 ) and (d 4 ). In all three cases one 
of the two inverse propagators generated by the WI will give rise to an expression similar 
to the second term on the RHS of ()2.3H) . i.e. a g-independent integral with a free Lorentz 
index, which cannot be saturated. These three terms are directly set to zero. The terms 



16 



stemming from the other inverse propagator read: 

q)Tf(k+£,q-£,-k-q), 
q)Tf(e + k-q,q-e,-k), 
q)Tf(k+t, -t-k). 

(2.34) 

and, as announced, can be directly identified with the corresponding terms on the RHS of 
(l2~m Therefore, 

<? (nj(9)| dl + nS(ff)U + nJ(?)L + n£fa)L) = o. (2.35) 

This concludes the non-perturbative proof of the special transversality property of IT^, (g): 
gluon and ghost loops are separately transverse, and the loops of different order do not mix 



1 u K b M\ d2 = *9f bet JJidk]id£]r;T D(k + £)A^(k)D(£- 
<ffi±{q)\, = Wf bci f f[dk][d£}T a ;: x D(k + £)A^(k)D(£- 



^nj(g)| d4 = igf bxi jj[dk}[d£}T a ^D(k + e)A^(k)D(£ 



D. Towards a new SD series 

As explained for the first time in the upshot the PT is to eventually trade the conven- 
tional SD series for another, written in terms of the new, gauge-independent building blocks. 
Then one could truncate this new series, by keeping only a few terms in a "dressed-loop" 
expansion, and maintain exact gauge-invariance, while at the same time accommodating 
non-perturbative effects. As mentioned in the Introduction, one of the most central issues in 
this context is how to convert the SD-series defining A^ into a dynamical equation, namely 
one that contains A^ v on both sides. For this to become possible, one must carry out inside 
the loops of the diagrams shown in the previous subsection the substitution A^ u — > A^ v . In 
this subsection we focus only on one particular aspect of this problem, namely the role that 
the BQIs might play in implementing the aforementioned substitution. 



The connection 39| between the PT and the Batalin-Vilkovisky quantization formalism 
55! ] has been instrumental for various of the recent developments in the PT. Specifically, 
using the formulation of the BFM within this latter formalism, one can derive non-trivial 
identities (BQI's), relating the BFM n-point functions to the corresponding conventional 
n-point functions in the covariant renormalizable gauges, to all orders in perturbation the- 



ory [44]. The relation between these two types of n-point functions is written in a closed 
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form by means of a set of auxiliary Green's functions involving anti-fields and background 
sources, introduced in the BFM formulation. These latter Green's functions are in turn 
related by means of a SD type of equation 39j to the conventional ghost Green's functions 



57 



appearing in the STI satisfied by the conventional all-order three-gluon vertex |56 , 

We will restrict our discussion to the case of the propagators and A M!/ . The relevant 
quantity appearing in the corresponding BQI is the following two-point function, to be 
denoted by A Qj g(g), defined as (we suppress color indices) 

KM = f [dk]H^D(k)A^(k + q) H; u (k + q, —k, -g), (2.36) 
where the elementary vertex H^l is 



H 



(o) 




-W9af3 



and H a/3 is given by 



H af )(p, r,q) 



H 



(o) 



a/3 



+ 




^2.37) 



with q + p + r = 0. JC u p is the conventional one-particle irreducible connected ghost-ghost 



gluon-gluon kernel appearing in the QCD skeleton expansion 
in the all-order Slavnov- Taylor identity satisfied by the conventional three-j 
and is related to the conventional gluon-ghost vertex ITg(p, r, q) by [l. Isfil I 

q a H af3 (p, r, q) = T p (p, r, q) . 



(2.38) 



Using the diagrammatic definition of H a p shown in Eq. ([2.37|h we may recast Eq. (j2.36J) in a 
dressed-loop expansion as follows, 




Due to the transversality of both li^q) and n^g), if we write 

Kp(q) = g a pG(q 2 ) + q a q p L(q 2 ), 



(2.39) 



(2.40) 
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it turns out that the fundamental all-order BQI between A(g) and A(g) involves only G(q 2 ), 



and is given by 



39 



44| 

A(q) = [1 + G(q)} 2 A(q). (2.41) 



It is then elementary to demonstrate that the full propagators (in the Feynman gauge) are 
related by 

A, v {q) = [1 + G{q)f A, v {q) - iG{q 2 ) [1 + G{q)\ ^ . (2.42) 

The process of replacing A^ v — > A^ v will therefore introduce the function G inside the 
loops; however, the theoretical and practical consequences of this operation are not clear to 
us at this point. Some of the various possibilities that one might envisage include: to study 
the dynamical equation for G {viz. Eq. (j2.36J) ) together with the SD of A, as a coupled 



system 



591 ]; attempt to reabsorb the G's into a redefinition of the vertices appearing in the 



diagrams, together with their corresponding SD equations; consider the diagrams containing 
G's as being of higher order in the dressed loop expansion, due to the additional explicit 
integration appearing in their definition, Eq. (j2.3(i|) . 

In the rest of this article we will adopt what appears to be the lowest order approximation 
in this context, setting A^ u = A^ u in the SD equation, and G = everywhere else. 



III. GENERAL CONSIDERATIONS FOR IR FINITE SOLUTIONS 

In this section we will briefly review some of the main issues involved when trying to 
obtain from the corresponding SD equation IR finite solutions for the gluon self-energy A. 
By "IR finite" we mean solutions for which A _1 (0) 7^ 0, which is equivalent to saying that 
A(0) is finite. 



A. Two necessary conditions 

Let us start by considering two necessary conditions for obtaining IR finite solutions. 

(i) From the one-loop dressed SD equation for the gluon self-energy (see Fig. |2J) it is 
clear that, since there can be at most two full gluon self-energies inside the diagrams, on 
dimensional grounds the value of A _1 (0) will in general be proportional to two types of 



19 



seagull-like contributions, 

To = J [dk] A(k) , 

Ti = / [rffc]A; 2 A 2 (A;). (3.1) 



Perturbatively, both % and 7i vanish by virtue of the dimensional regularization result 

!^lln^ 2 )=0, N= 0,1,2,... (3.2) 

which guarantees the masslessness of the gluon to all orders in perturbation theory. In order 
to permit IR finite solutions one must assume that seagull-like contributions, such as those 
shown in (J3.1|) . do not vanish non-perturbatively. Of course, once assumed non- vanishing, 
both T and T\ are quadratically divergent, and, in order to make sense out of them, a 
suitable regularization must be employed. 

(ii) The form of the full three-gluon vertex is instrumental for the generation of IR 
finite solutions. Specifically, as is well known from the classic papers on dynamical mass 
generation 4^|, in order to obtain A _1 (0) 7^ one must introduce massless poles in the full 
three-gluon vertex F„ Q/ 3, appearing in the expression for U^\q), Eq. (j2.12j ) [go] ]. Notice in 
particular that, whereas after allowing for non- vanishing seagull contributions the inclusion 
of graph (a 2 ) is essential for the transversality of IT^, its presence does not lead to A _1 (0) 7^ 
0. Thus, if the full three-gluon vertex IT satisfies the WI of (12.22}) . but does not contain 
poles, then the seagull contribution T 7^ of graph (a 2 ) will cancel exactly against analogous 
contributions contained in graph (ai), forcing A _1 (0) = 0. Put in a different way, the non- 
vanishing seagull contribution that will determine the value of A _1 (0) is not the one coming 
from graph (a 2 ). This is even more evident in the non-linear treatment, where the term Ti 
makes its appearance; clearly, such a term cannot be possibly obtained from (a 2 ). 

To appreciate the delicate interplay between the points mentioned above, let us consider 
the ghostless, one-loop dressed version of the SD equation, by keeping only the two graphs 
of group (a). The SD equation has the general form, 



with 



a- V)P^0?) = g 2 PM + * + rrj, 2) ] , (3.3) 

n^O?) = \c A g 2 J[dk]f^K aa \k)f ua ,^\k + q) ) 

n£ 2) = -C A g 2 g,u%. (3.4) 
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Let us write zTljL* (g) in the general form 



fifc\q) = q 2 A(q 2 )g^ u + B(q 2 )q^q v , (3.5) 



where A(q 2 ) and B(q 2 ) are arbitrary dimensionless functions, whose precise expressions 
depend on the details of the Jr va 'p' employed. 



The transversality of i 



ii£ ) ( g ) + n<E ,) 



implies immediately the condition 
q 2 [A{q 2 ) + B{q 2 )] = iC A g 2 %, (3.6) 
and thus the sum of the two graphs reads 

i [u^\q) + n£>] = -q 2 B{q 2 )P, u {q) • (3-7) 

Clearly, 

A- x (0) = \im(-q 2 B(q 2 )) = lim (q 2 A(q 2 )) - iC A g 2 % . (3.8) 

Interestingly enough, once the transversality of n M „ has been enforced, the behavior 
of A _1 (g 2 ) is determined solely by B(q 2 ). In particular, the value of A -1 (0) is given 
by lim (—q 2 B(q 2 )). Evidently, if B(q 2 ) does not contain (1/q 2 ) terms, one has that 

lim (—q 2 B(q 2 )) = 0, and therefore A _1 (0) = 0, despite the fact that T has been assumed 
to be non-vanishing. Actually, as we will see in the context of the linearized SD equation 
that we will study in the next section, if does not contain massless poles, it is precisely 
this latter situation that is realized, by virtue of an identity relating the various integrals in- 
volved. On the other hand, if B(q 2 ) contains (1/q 2 ) terms, lim (—q 2 B(q 2 )) ^ 0, allowing for 
A _1 (0) 0. In fact, for physically acceptable solutions, one must demand that A _1 (0) > 0, 
which imposes further restrictions on the possible forms of IT. Of course, this is not to say 
that the presence of poles in the vertex is sufficient for obtaining IR-finite solutions, because 
B(q 2 ) may end up being non-singular due to accidental algebraic cancellations. As we will 
see in a concrete example in the next section, the net contributions of pole terms originating 
from different Lorentz structures within the same vertex may lead to an equation that does 
not generate mass. 

The quantities A(q 2 ) and B(q 2 ) appearing in (|3.5)1 . will be functionals of the unknown 
quantity A, as dictated by the SD equation; their specific form will depend on the details 
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of the vertex Ansatz chosen. In general, their value at q 2 = will be linear combinations of 
the two terms defined in (|3.1|) . namely 

A-^O) = g 2 (a % + aiTi) , (3.9) 

where the values of the numerical coefficients and a\ depend on the details of the problem. 
Then, the corresponding SD equation will read schematically 

AT V) = A- x (0) + g 2 q 2 J [dk] F{A; q, k} , (3.10) 

where the functional J^{A; q, k} is regular at q 2 = 0. Then, after renormalizing, we will be 
looking for solutions of (J3.10)) that simultaneously: (a) reproduce correctly the asymptotic 
behavior for A(q 2 ) predicted by the RG, viz. Eqs. (|2.9|) and (12. (b) are finite at q 2 = 0; 
(c) satisfy the (appropriately regulated) constraint of ()3.9|) [6lJ. Turns out that points (a), 
(b) , and (c) are deeply intertwined, in a way that we will sketch below, and further elaborate 
upon in Sec. IIVI 



B. RG behavior and the SD equation 

The asymptotic behavior that A(g 2 ) must satisfy in the deep UV is given by Eqs. (J2.9|) 
and (j2.11|) . In practice, however, it is highly non-trivial to obtain, from the corresponding 
SD equation, solutions displaying this asymptotic behavior. This difficulty is intimately 
connected to the approximations used for the (all-order) vertex IT. Of course, within a full 
SD equation treatment, IT satisfies its own non-linear integral equation, which determines 
its structure. One must deal then with a very complex system of coupled integral equations 
involving A, IT, and several many-particle kernels. The usual way to reduce the difficulty of 
this problem is to resort to the gauge-technique, namely express F as a functional of A, in 
such a way as to satisfy (by construction) the first WI of Eq. ()2.22|) exactly. This procedure 
fixes the "longitudinal" part of the vertex, but leaves its "transverse" (identically conserved) 
part undetermined. This ambiguity, in turn, leads to the mishandling of overlapping diver- 
gences, which manifests itself in the fact that (i) one cannot renormalize multiplicatively, 
but only subtractively, and (ii) the RG-behavior of the solutions is distorted. In particular, 
one obtains solutions which asymptotically behave like 

A-\q 2 )=q 2 



q 

1 + wbg z In I — 
ix 2 



(3.11) 
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with w > 1. These solutions reproduce upon expansion the expected (one- loop) perturbative 
result, but non-perturbatively they miss the correct RG behavior of ([2.9)1 . 

The first-principle remedy of the situation would require the full treatment of the con- 
served part of the vertex, in a way similar to that followed in 62] for the vertex appearing in 
the electron and quark gap equations. Unfortunately, extending their method to the case of 
the three-gluon vertex is technically very involved, and is at the moment beyond our powers. 
Instead, we propose to model the RG behavior according to the simple prescription put forth 
in 11, 63[. The basic observation is that the correct RGI may be restored if every A(z) 
appearing inside ^{A; q, k} were to be multiplied ("by hand") by a factor (see Eq. (|2.11j )) 



9 



2 



l + bg 2 \n(z/^) = ^--. (3.12) 

9 \ z ) 

Thus, one is effectively switching from (|3.1Uj) to the corresponding "RG-improved" equation 

A- V) = A- x (0) + g 2 q 2 J [dk] F{A(z) (g 2 /g 2 (z)) ;q,k}, (3.13) 

with z = k 2 or z = (k + q) 2 ; equivalently, in terms of the manifestly RG-invariant quantities 
d(q 2 ) and g 2 (q 2 ), 



d~\q 2 ) = d~\0) + q 2 J[dk] F{d(z)/g 2 (z);q,k} (3.14) 

and 

<H(0) = a Q T Q + a{T x , (3.15) 
where (in Euclidean space), by virtue of ()2.10|) 

d(k 2 ) 



To = J [dk]j 



g 2 (k 2 ) ' 
k 2 d 2 {k 2 ^ 
¥(k 2 



C. Regularization of seagull-like terms 



Returning the issue of the regulation of the constraint ()3.15j) . notice that the integrals 
in (|3.1|) differ from those of (|3.16J) in an important way. Roughly speaking, the presence 
of the RG logarithms in their numerators, (contained in g 2 (k 2 ) and g 4 (k 2 ), respectively) 
compensates the logarithms contained in d(k 2 ) and d 2 (k 2 ), allowing one to regularize them 
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simply by subtracting a unique (vanishing) integral, that of Eq. (j3.2j) for N = 0, provided 
the solutions satisfy certain generic conditions. 

To study this in detail, let us for the moment concentrate on (potential) solutions of the 
SD equation that are qualitatively of the general form 

d(q 2 )= 2 f p( C } 2r (3.17) 
q A + m z (q z ) 

where 



gl P (q 2 



bin 



q 2 + f(q 2 ,m 2 (q 2 )) 



(3.18) 



A 2 

The function m 2 (q 2 ) may be interpreted as a momentum dependent "mass" with the property 
(to be imposed self-consistent ly) that m 2 (0) > 0. In addition, we expect that m 2 (q 2 ) is a 
monotonically decreasing function of q 2 , with m 2 (q 2 ) — > as q 2 /A 2 — ► oo. The quantity 
^np(<? 2 ) represents a non-perturbative version of the RG-invariant effective charge of QCD, 
going over to g 2 (q 2 ) in the deep UV. The (dimensionfull) function f(q 2 , m 2 (q 2 )) is expected to 
be such that 'g 2 iV ,{q 2 ) will be a monotonically decreasing function of q 2 , with /(0, m 2 (0)) > A 2 . 
The dimensionality of / is to be saturated by m 2 (q 2 ); thus if one were to set m 2 (q 2 ) = 
then one should have f{q 2 ,0) = 0. The presence of a f(q 2 ,m 2 (q 2 )) with such properties 
in the logarithm of ^ p (g 2 ) eliminates the Landau pole, and leads in the deep IR to the 
characteristic property of "freezing" . For the analysis that follows, note also that 

1 1 ( f(q 2 ,m 2 (q 2 ))\ 

bin 1 + - { - - V - 11 . (3.19) 



9np(? 2 ) g 2 (q 2 ) V q 2 

It turns out that, for the proposed regularization to work, both m 2 (k 2 ) and / (q 2 , m 2 {q 2 )) 
must drop "sufficiently fast" in the deep UV. 

To understand this point, we substitute into T of Eq. ffi.lfij) a solution of the form ()3.17j) . 



^o= l[dk] m , ^> h2 ^ (3-20) 



[k 2 + m 2 (k 2 )]g 2 (k 
and consider T™ s , obtained after subtracting j[dk]/k 2 = from T , 



r reg - !\dk] ( 9lp{e) - -) 

° ~ J \[k 2 + m 2 {k 2 )]g 2 {k 2 ) k 2 J 



Let us examine the two integrals on the RHS separately. If m 2 (k 2 ) behaved asymptotically as 
ln~ a (k 2 ), with the anomalous mass-dimension a > 1, then the first integral would converge, 
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by virtue of the elementary result 

dz 1 



/ 



(3.22) 



z(\nz) 1+ ~' 7(1hz)t' 
Whether such a behavior of m 2 (k 2 ) is realized or not must be verified directly from the 
corres pon ding SD equation. For example, this was indeed the case for the equations studied 
in [ill Ig.sI . (with a = 12/11), and we will observe it again in the next section. In fact, a faster 
asymptotic behavior of the form \n a (k 2 )/k 2 may be obtained from non-linear versions of the 
SD equation The second integral will converge as well, provided that f(k 2 ,m 2 (k 2 )) 

drops asymptotically at least as fast as ln~ c (/c 2 ), with c > 0. If, for example, / = pm 2 {k 2 ) 
(with a > 1, for the first integral to converge), then the convergence condition for the second 
integral is automatically fulfilled. Notice that perturbatively T™ 5 vanishes; this is because 
m 2 (k 2 ) = to all orders, and therefore, since in that case also / = 0, both integrals on the 
RHS of (EEZQ) vanish. 

Assuming that m 2 (k 2 ) and f(k 2 ,m 2 (k 2 )) behave as described above, then it is straight- 
forward to verify that the difference T\ — T is automatically finite. Indeed, applying the 
elementary identity k 2 = [k 2 + m 2 (k 2 )] — m 2 (k 2 ) in the numerator of T l5 we arrive at 

* - T. = - / [*, - 6 /[*, fg^) In (x + , ,3,3, 

where both integrals on the RHS converge, without any additional assumptions. 
Thus, the RHS of ()3.15|) can be written as 

d-\0) = (a + ai)T + a l {T 1 - T ) . (3.24) 

Clearly, if we happened to have that the RHS of (|3.24j) would be automatically 

convergent, without further need of regularization. If a\ ^ — ao, we will replace on the RHS 
of ()3.24|) T ) by T™ 9 , arriving at the regularized version of (|3.15|) 

Cg(0) = (ao + ai)T rcg + a 1 (T 1 - T ) . (3.25) 

It is important to emphasize that the form of the solutions assumed in (|3.17|) is meant 
to quantify the UV behavior necessary for the proposed regularization to work, but does 
not restrict their deep IR behavior. Specifically, let us assume that a solution of the type 
()3.17|) . to be denoted by d c (q 2 ), satisfies the necessary asymptotic conditions, and consider 
a new function, d-^(q 2 ) = dc(q 2 ) + h(q 2 ), where h(q 2 ) vanishes faster than d c (q 2 ) in the 
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UV, but is not restricted in the IR. Then, if d-^(q 2 ) is inserted into (J3.21j) . the corresponding 
integrals are still finite. As we will see in Sec. this situation does in fact occur: when 
solving the corresponding SD equation, in addition to the "canonical" solution of the type 
f!3.17|) . we find solutions that in the UV go over to f)3 . 1 7j) . but in the deep IR display a much 
sharper increase. The point is that these new solutions can be regulated following the same 
procedure outlined here. 

Several comments are now in order. 

(i) The above method for regulating the seagull-like terms relies solely on the integration 
rules of the only known gauge-invariant regularization scheme, namely dimensional regu- 
larization, together with the requirement of an appropriate momentum dependence for the 
dynamical mass. In that sense it is conceptually rather economical, evoking a minimum 
amount of additional theoretical input. 

(ii) The implementation of the proposed regularization hinges crucially on the requirement 
that, within the given truncation scheme, the RGI behavior can be encoded faithfully into 
the SD equation, and the logarithmic terms are correctly accounted for. In particular, 
the compensation (in the UV) of the RG logarithm contained in d(k 2 ) by the logarithm 
of g~ 2 (/c 2 ), is essential for the consistency of the ensuing regularization of (j3.16j) . In the 
absence of the compensating logarithm one would have to subtract instead a term 1/k 2 In k 2 
in order to achieve UV convergence. The rules of dimensional regularization allow such a 
possibility; by virtue of the more general result 3 

[dk](k 2 ) a = 0, (3.26) 



valid for any value of a, together with the elementary identity 

ml dte ~ ttR '^ k2 ^ 2 r tb92 (3-27) 



l + bg 2 \n{k 2 /fi 2 " 



(valid for R > 0), one may set 

[dk] 



dt e~* (/i 2 )^ 2 [dk] (k 2 )- {1+tb92) = 0. (3.28) 



k 2 [l + bg 2 \n(k 2 /p?)\ 

Subracting such a term would eventually regulate the initial integral in the UV. The problem, 
however, is in the IR: the logarithm in the denominator of the regulated integral would give 
rise to the very pathology one has set out to cure in the first place, namely the Landau pole. 
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(iii) Of course, d~ 1 (0) must be positive definite (in Euclidean space). The regularization 
possibility offered by (|3.21|) was in fact appreciated in [llj], but was not pursued further, 
on the grounds of furnishing the "wrong" sign for d -1 (0). In the context of that work this 
was indeed so, because the three-gluon vertex used (the analogue of IT) was completely 
fixed, being the tree-level vertex derived from the Lagrangian of the massive gauge invariant 
Yang-Mills model. After connecting T to the finite vacuum expectation value of a composite 
scalar field creating + glueball states, the replacement 

J [dk] A(k 2 ) — > A-\0) J [dk] A 2 (k 2 ) (3.29) 

was used instead. Assuming that the correct RG-behavior is captured by the corresponding 
SD equation, then the integral on the RHS of ()3.29|) converges due to the extra logarithms 
in the denominator. 

In the case we consider here, the sign situation is more involved. The two integrals on the 
RHS of Eqs. (j3.21|) are positive definite (assuming that m 2 (k 2 ) > and f(k 2 ,m 2 (k 2 )) > 
in the full range of momenta); thus, the sign of Tq GS is fixed. On the other hand, the sign 
of the RHS of ()3.23j) is not definite, since the In z appearing in the numerator of the second 
integral (contained in g~ 2 (A; 2 )) becomes negative at z < 1. In addition, and perhaps more 
importantly, the signs of ao and a± are not a-priori known either. This is so because the 
expression for <i _1 (0) is determined not by the (fixed and known) sign of the seagull graph, 
but from a delicate combination of the coefficients of the pole terms contained in the full 
three-gluon vertex T. In our opinion the issue of the sign should be settled within the strict 
confines of QCD and dimensional regularization. Thus, if a QCD-derived (approximate, to 
be sure) form for IT were to yield (after applying the proposed regularization) a negative 
sign for d _1 (0), then one should be inclined to conclude that the dynamical mass generation 
is not realized, at least not in the context of the specific truncation scheme. 



IV. THE LINEARIZED SD EQUATION 

In this section we will study in detail a linearized version of the SD equation obtained in 
the one-loop dressed approximation, omitting the ghosts. The resulting equation is a close 



variant of the one presented in 
address further points of interest 



but displays several distinct features, allowing us to 
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A. Linearizing the SD equation 



We start by considering the two diagrams of group (a), shown in Fig. El Since we are 
working in the Feynman gauge of the renormalizable gauges, instead of the axial gauges 
used in jllj], the general form of A^ is that of Eq. ([2.4)1 . In order to be able to use the 
simplified Ansatz for the vertex given below in conjunction with the Lehmann representation, 
it is necessary to drop the longitudinal parts of A^ inside the integrals, using A fJiU (k) = 
—ig^ w A(k). As we will see in a moment, omitting these terms does not interfere with the 
transversality of the external A fiu (q); in a way it is like considering scalar QED, with massive 
scalars inside the vacuum polarization loop, yielding a transverse photon self-energy. 

After dropping the longitudinal parts inside the loops, we obtain 

ff^(g) = \c A g 2 ^ J[dk]ffA(k)f u ^A(k + q) - 2dg^J[dk] A(k)J , (4.1) 

with 

?W; = (2k + q)^g a p - 2q a g^ + 2q /3 g IMX , (4.2) 

and 



vaj3 



A-\k + q)-A-\k) g Ql3 . (4.3) 



Then it is straightforward to check that IL^g) is transverse. 



As in llj, in order to linearize the SD equation of (j4.1|) . we resort to the Lehmann 
representation for the gluon propagator, setting 



P(A 2 ^ 
q 2 — X 2 + it 



This way of writing A(q 2 ) allows for a relatively simple gauge-technique Ansatz for T^ a g, 
which linearizes the resulting SD equation. In particular, on the RHS of the first integral in 
(|4.1|) one sets 

where T^ a g must be such as to satisfy the tree-level WI 

q"? ^ = [(k + q) 2 - k 2 } g aB = [(k + qf - A 2 ] g aB - (k 2 - X 2 )g aj3 . (4.6) 

Then it is straightforward to show by contracting both sides of ()4.5|) with q u , and employing 
fO)J) and (Q, that satisfies the all-order WI of Ea. ljOjl . Of course, choosing V ua3 — 
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T ua p solves the WI, but as we will see in detail in what follows, due to the absence of pole 
terms it does not allow for mass generation, in accordance with the discussion in the previous 
section. Instead we propose the following form 



The essential feature of this Ansatz is that, due to the inclusion of the 1/q 2 pole term, it can 
give rise to IR finite solutions. Note that the additional terms have the correct properties 
under Bose symmetry with respect to the two quantum legs. For our purposes the constants 
Ci, c 2 , and c 3 are treated as arbitrary parameters, offering the possibility of quantitatively 
examining the sensitivity of the solutions on the specific details of the form of the vertex. 
Of course, in reality their value will be determined by the dynamics of the corresponding 
SD equation satisfied by the full vertex, a problem which is beyond our powers at present. 
Let us next define the quantities 



Substituting in ()3.4|) the expression for F obtained from the combination of ()4.5|) and (j4.7|) . 
after some straightforward algebra we obtain 




(4.7) 




(4.8) 




(4.9) 



with 



F(q 2 ,\ 2 ) = (7d-8)q 2 B(q 2 ,\ 2 ) + 2d\2\ 2 B(q 2 ,\ 2 )-(d-2)T(X 2 ) 



+da - q 2 B(q 2 , A 2 ) + A\ 2 B(q 2 , A 2 ) + 2T(A 2 ) 



+4(rf - l)c 2 \ 2 B(q 2 , A 2 ) + T(A 2 ) + 4(d - 1) c 3 q 2 B(q 2 , A 2 ) . (4.10) 



In order to study whether A x (0) 7^ 0, we must determine the value of F(0, A 2 ). To that 
end, note the crucial identity 



2A 2 £(0,A 2 ) = (rf-2)T(A 2 ), 



(4.11) 
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which may be easily proved following the standard integration rules of dimensional reg- 
ularization. Thus, if we were to eliminate the pole terms in the vertex (|4.7jl by setting 
Ci = c 2 = 0, then, by virtue of (|4.11j) . we have F(0, A 2 ) = 0. Evidently, the terms propor- 
tional to Ci and c 2 in (J4.1(Jj) are non- vanishing at q 2 = 0, even after the application of (J4.ilj) . 
thus yielding A _1 (0) ^ 0. To see this in detail, let us define 

B(q 2 ,X 2 ) = B(q 2 ,X 2 )-B(0,X 2 ), (4.12) 

and then replace everywhere in ()4.10j) B(q 2 ,X 2 ) = B(q 2 ,X 2 ) + 5(0, A 2 ), using ()4.11j) to 
eliminate X 2 B(0, A 2 ) in favor of T(A 2 ). This leads to 

F(g 2 ,A 2 ) = g 2 |(7d-8)-rfc 1 + 4(rf-l)c 3 |^(0,A 2 ) + 5(g 2 ,A 2 

+ 4rf(l + ci) + (d - l)c 2 X 2 B(q 2 , A 2 ) + 2d(d- 1) (ci + c 2 ) T(A 2 ) , (4.13) 
from which follows immediately that 

A-\0) =d( Cl + c 2 )C A g 2 i J[dk]A(k 2 ). (4.14) 

Thus, the value of A _1 (0) is determined solely from the singular part of graph (cxi), in 
agreement with the general discussion of the previous section. In particular, the seagull term 
corresponding to graph (a 2 ), due to the aforementioned cancellation, imposed by (|4.11|) . does 
not enter in the expression for A _1 (0). This is of course not to say that (a 2 ) is irrelevant; 
on the contrary, as we have seen, the role of (a 2 ) is crucial in enforcing transversality, 
which eventually allows one to arrive at Eq. (j3.8|) and Eq. (j4.14J) . Notice also that the factor 
determining the value of A~ 1 (0) is the sum c\ + c 2 ; thus, one could envisage the possibility 
of contributions from pole terms pertaining to different Lorentz structures canceling against 
each other, or yielding the wrong sign for A _1 (0) . 

B. Further algebraic manipulations 

We will now further manipulate Eq. ()4.13|) . The term proportional to B(0, A 2 ) on the RHS 
of ()4.13p diverges, and is to be absorbed into the wave-function renormalization constant, 
soon to be introduced. Using the expression for 5(0, A 2 ) obtained directly from ()4.8|) . we 
can write it in the alternative form 
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where we have assumed that the order of integration may be changed. In addition, we will 
use the elementary result 



BUr.X-^^j ^Inh+i^A) (1.16) 



together with the following identities 



dX 2 p(X 2 ) l 1 dx\n(l + q2x{X x2 1} ) = ^^^l-^y^A^) 
1 / q 2 x(x-l)\ f q2/4 , / Az^ 1/2 



o 



dX 2 p(X 2 )X 2 j dxlnl 1 + - ~ v ~ 2 = y dzz\^l-^j A(z), (4.17) 

which allow us to rewrite the RHS manifestly in terms of the unknown function A. At this 

point it is obvious that the perturbative result, which must be proportional to (7d — 8) [see 

n^(g) in flUBD ], will be distorted by the presence of the other two terms inside the curly 

brackets on the RHS of (j4.13|) . To avoid this we will use the freedom in choosing the value 

of c 3 , and fix it such that 4(d— l)c 3 = dc\. After this, using the results given above, setting 

d = 4 everywhere except in the measure, and defining 

~ 10C A 6( Cl + c 2 ) 4 + 4 Cl + 3c 2 

6 "48^' 7 " 5 ' (418) 

we arrive at the integral equation 



, q2/A ( Az\ l/2 - 
+ 1 bg 2 I dzz[\- — \ A(z) +A- 1 (0), (4.19) 

with 



o V <? 2 



A-i(o) = ^ ld d kA(k 2 ). ( 1.20) 



Next consider the Euclidean version of (|4.19|) : to that end we set q 2 = —q 2 , with q 2 > 
the positive square of a Euclidean four-vector, define the Euclidean propagator as 

A E (q 2 E ) = -A(-q 2 E ) : (4.21) 

and the integration measure [dk] = i[dk] E = id k E /(2it) . To avoid notational clutter, we 
will suppress the subscript "E" everywhere except in the d d k measure. Then we have 



+ lbg 2 dzzil- — ) A{z)+A- 1 {0), (4.22) 



° 2/ \ (. 4 £ \ 1/2 
V Q 2 
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where from now on q 



2 stands for the (positive) square of a Euclidean vector, and 



A" 1 



(0) 




/ 



d d k E A(k 2 ). 



(4.23) 



C. Renormalization 

In order to renormalize the equation, first we define the bare and renormalized quantities 
as follows: 



by virtue of the Abelian-type Wis satisfied. Then, it is straightforward to verify that the 
net effect of renormalizing (j4.1J) . or subsequently f)4.19|) . is to simply multiply its RHS by 
Za and replace all bare quantities by renormalized ones. 

However, as usually happens at this level of approximation, where the overlapping diver- 
gences are not properly accounted for, due to the ambiguities in the longitudinal parts of 
one is forced to renormalize subtractively instead of multiplicatively. This amounts to 
interpreting Za as an infinite constant that renders the product 





(4.25) 



finite, and setting Za = 1 in all other terms. This procedure leads to 




(4.26) 



The renormalization constant K is to be fixed by the condition 



a-V) = /z 2 , 



(4.27) 



with /i 2 a Euclidean momentum, satisfying fi 2 3> A 2 , 



yielding 




(4.28) 
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D. Renormalization-group analysis 



We next study the UV behavior predicted by the integral equation ()4.26|) for A(q 2 ). To 
begin with, notice that the perturbative one-loop result may be recovered by replacing on 
the RHS of (|4.26j) A(z) — > 1/z, the tree-level value; then, the second term vanishes, and 
after setting (1 — Az/q 2 ) 1 ^ 2 — > 1 and (1 — Az/^i 2 ) 1 ^ 2 — > 1 in the first term (curly brackets), 
we obtain A _1 (g 2 )| pert = q 2 {l + bg 2 ln(g 2 / /i 2 ) J . However, if one were to solve this equation 
non-perturbatively, one would discover that, even though the perturbative result is correctly 
recovered, A(q 2 ) does not display the expected RG behavior, e.g. that of Eq. tj2.9j) . with 
b^b. 

The fact that this behavior is not captured by (|4.26|) can be easily seen by writing down 
the simplified version of that equation, 



A-V) = q 2 \l + bg 2 / dzA(z) , (4.29) 




valid in the deep UV, and converting it into an equivalent differential equation. Setting 
A(g 2 ) = G(q 2 )/q 2 , and x = q 2 , we obtain 

i_m = _ lg ,9M, (4.30) 

which leads to 



A- 1 (q 2 ) = q 



l+2bg 2 In 



f' 2 



(4.31) 



Upon expansion this expression yields A l (q 2 )\ per t correctly, but differs from (J2.9J) . The 
fundamental reason for this discrepancy can be essentially traced back to having carried out 
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63j , thus distorting the RG 



the renormalization subractively instead of multiplicatively 
structure of the equation. 

In order to restore the correct RG behavior at the level of (|4.2fij) . we will use the procedure 
explained in Sec. \\\\\ substituting in the integrands on the RHS of (|4.2fij) 

AW^^f. (4,2) 
Then (|4.26J) may be rewritten in terms of two RG-invariant quantities, d(q 2 ) and g 2 (g 2 ), as 
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follows 



I _ r q2/i ( Az\ 1/2 
d-\tf) = q 2 lK' + bj dz fl-^J 



Az\' A d(z) 



with 



and 



~ f q2/4 ( Az\ 1/2 d(z) 

+ 76 / dzz (l- — ) +d~\0), (4.33) 

Jo \ TJ 9 (z) 

, 1 r T 2/4 , ( z \ ( Az\ 1/2 d(z) , , 



7T Z 



d _i (0) = z d a k B -±—^. (4.35) 



S 2 (P) 



It is easy to see now that Eg. ()4.33j) yields the correct UV behavior for d(q 2 ), given in (|2.1()jl . 
For example, converting ()4.33|) into a differential equation, and setting d(q 2 ) = F(q 2 )/q 2 , 
the equivalent of ()4.30|) now reads 

^M = -6^, (4.36) 
dx x g (x) 

whose solution is F(x) = g 2 (x), as announced. 



E. Asymptotic behavior of m 2 (q 2 ) 

As has been discussed in Sec. IHIl the ability to regularize condition ()4.35|) following the 
properties of dimensional regularization depends crucially on the asymptotic behavior of 
m 2 (q 2 ). It is therefore essential to determine the asymptotic behavior that Eq. (J4.33|) yields 
for m 2 (q 2 ) in the deep UV. To that end, substitute (j3.17j) on both sides of (j4.33jl . and 
consider the limit where q 2 is large. The integral proportional to 7 may be written as 

fll iz U _ <£) V! __ = t _ f'\ z U - *!) V2 . (4.37) 

\ q 2 J z + m 2 {z) 6 J Q \ q 2 J z + m 2 (z) 



Then setting (1 — Az/q 2 ) 1 ^ 2 — > 1, and dropping terms that do not grow logarithmically with 
q 2 , f|4.33|) reduces to 

f q2/A dz f q2/4 m 2 (z) 

[q 2 + m 2 (q 2 )] In (q 2 /A 2 ) = q 2 / — ™— - 7 / cfe-£-W . (4.38) 

Jo z + m 2 (z) Jo z + m 2 (z) 

Let us then separate the terms that go like q 2 \nq 2 and m 2 (q 2 ) In q 2 ; obviously the first 

integral on the RHS compensates the q 2 \nq 2 on the LHS. Then setting 

m 2 (q 2 ) ~m 2 \n- a (q 2 /A 2 ) , (4.39) 
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we find that the two sides of the equation can be made equal if 

a = 1 + 7. (4.40) 

Thus, provided that 7 > 0, the momentum dependence of the mass term in the deep UV is 
of the type needed for the regularization of Eq. (j3.21|) to go through. 



V. NUMERICAL ANALYSIS 

The equation we will solve is ()4.33|) with the renormalization constant K' of ()4.34|) ex- 
plicitly incorporated, e.g. 




az\ 1/2 d{z)_ ^_ r 2/ * c 

q 2 J g 2 (z) [i 2 J 



Az 


\ 1/2 d(z) 


'J 2 


) 9 2 (z) 


Az\ 


1/2 d{z) N 


J 2 ) 


9 2 {z) , 



(5.1) 



with d 1 (0) given by ()4.35|) . eventually to be replaced by its regularized expression, according 
to (|3.21j) . In particular, if the numerical solution obtained for d(q 2 ) has the general form of 
Eq.(jSHZ!), d-\0) will be given by 



d-\0) = la 



m 2 (z) ~ f°° J 7 . . , / f(z,m 2 (z)) 



dz r— + b / dzz d(z) In 1 + 



(5.2) 



Uo z + m 2 (z) J " ' " V z 
If the solutions deviate in the IR from Eq. (j3.17|) . an analogous expression can be obtained, 
(see below) in accordance with the discussion following Eq. (j3.25j) . 

Equations (j5.1j) and (|5.2jl form a system of equations that must be solved simultaneously. 
The role of (|5.1|) is to furnish possible solutions for d(q 2 ), while f!5.2|) constrains them or 
the value of the parameters involved. Roughly speaking, the strategy for solving the system 
is the following. Note first that (|5.1|) contains no information for the value of <i -1 (0); at 
q 2 = 0, one obtains an identity. Therefore, one chooses an arbitrary value for d _1 (0), 
plugs it in on the RHS of (|5.1j) . and proceeds to solve the integral equation. Given such 
a solution, one must then substitute it in ()5.2j) and calculate the value of <i _1 (0) obtained 
after the integration, and compare it with the value of <i _1 (0) chosen at the beginning, the 
objective being to reach coincidence between the two values. Assuming that the value of 
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is considered to be fixed, one must then repeat this procedure, varying the initial value of 
<i _1 (0), until agreement has been reached. If instead one is free to choose the value of a, 
then for a given initial value of c? _1 (0) one varies a until compliance has been achieved. In 
this article we will follow the latter philosophy, treating a as an adjustable parameter. 

Specifically, for each value of c? _1 (0) chosen, we should vary 7 and a, in order to scan the 
two-parameter space of solutions. In practice, we choose to reduce the number of parameters 
down to one, namely a, since from the ensuing numerical analysis it becomes clear the 
dependence of the solution shows a very mild dependence on 7. Therefore, we rewrite 7, 
given by Eq. (j4.18|) . as 

7-f* + !-f, (8.3) 

and then we set C2 = 0, in the order to keep a as the only free parameter. The small 
differences produced when C2 is non-vanishing will be commented later on. 

For the numerical treatment we define a logarithmic grid for the variables q 2 and z\ this 
improves the accuracy of the algorithm in the small q 2 region, since the size of the steps 
is made smaller for IR momenta. We split the grid into two region: [0, /x 2 ] and (/^ 2 ,Auv]- 
Such splitting in needed for imposing on d(q 2 ) the renormalization condition (given by 
Eq. (j4.27|) ) at a perturbative scale jj 2 . Typically, we chose /1 2 = M| = (91. 18) 2 GeV 2 and 
Ajjv = 10 6 GeV 2 . Furthermore, one has to specify the coupling g 2 = g 2 (/U 2 ) entering in 
Eq. (j5.1|) : its value is obtained from Eq. ()2.1H) . where we have properly replaced b — > b, and 
used as input a value of A = 300 MeV for the QCD mass scale. Then, we solve the integral 
equation iteratively, starting out with an initial trial function, and using as a convergence 
criterion that the relative difference between the input and the output should be smaller 
than 10 ~ 8 . 

The general trend displayed by the solutions is that the characteristic IR plateau ( "freez- 
ing") becomes increasingly narrower as one increases the value of d(0). We have found two 
types of solutions, depending on the initial value chosen for d _1 (0): (i) For values of <i" 1 (0) 
within the ranee [0.01 GeV 2 , 0.07 GeV 2 ] the solutions can be perfectly fitted by Eq. (jHHZD ; 
we will refer to them as "canonical", (ii) For d _1 (0) < 0.01 GeV 2 the solutions can be 
fitted by Eq. ()3.17|) until relatively low values of q 2 , but deviate significantly in the deep 
IR, where they display a sharp rise and a rather narrow plateau; we will call such solutions 
"non-canonical" . 
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A. Canonical solutions 



In Fig. we show a typical case of a canonical solution corresponding to the initial choice 
d -1 (0) = 0.04 GeV 2 . As can be observed from the plot, d(q 2 ) is essentially a constant, 
determined by d(0), until the neighbourhood of q 2 = 0.01 GeV 2 ; then, the curve bends 
downward in order to match with the perturbative scaling behavior at a scale of few GeV 2 . 
All such solutions may be fitted very accurately by means of the d(q 2 ) of Ea. (l3~T7l) . where 
the functional form of g 2 p(<? 2 ) is that of Eq. ()3.18|) . with the function f{q 2 ,m 2 (q 2 )) fixed as 

m i {q 2 ) 



f(q ,m (q )) = pim (q ) + p 2 



and the dynamical mass has the form 



with exponent a 



m 2 (q 2 ) 
1+7- 



m n 



In 



q 2 + m 2 



v q 2 + pim 2 > " 1 
A 2 



In 



Pi rap 
A 2 



(5.4) 



(5.5) 
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FIG. 6: The black dots represent the numerical solution for d{q 2 ) obtained for the choice d x (0) = 
0.04 GeV 2 , a = 0.212 and 7 = 0.941. The continuous red line is the fit of Eq. (pn7|) . setting 
p 2 = -3.208 and = 0.45 GeV 2 . These curves are obtained fixing d _1 (0) = 0.04 GeV 2 , 
a = 0.212 and 7 = 0.941. The dashed blue line is the one-loop perturbative behavior. 

The dynamical mass, m 2 (q 2 ), and the running charge, a(q 2 ) = ^ 2 p (g 2 )/47r, corresponding 
to the numerical solution presented on the Fig. El are shown in Figs. [7| and |H1 respectively. 
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FIG. 7: The dynamical mass, m 2 (q 2 ), corresponding to the solution of Fig. El 
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FIG. 8: The running charge, a(q 2 ) = g 2 p (q 2 )/ATr corresponding to the solution of Fig. El 



In fact, in all cases studied, pi was fixed to the value p\ — 4, which was found to minimize 
the x 2 °f t ne corresponding fits. Therefore, the unique free parameter is P2, since the value 
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of wig can be written in terms of p 2 , simply by setting q 2 = 0, in Eq. (j3.17j) . i.e., 

d~\0) = bml\n(l^^ , (5.6) 

where 

f(0,m 2 Q ) = (4 + p 2 )m 2 . (5.7) 

Thus, the value of the IR fixed point of the running coupling is determined by the value 
assumed by /(0, m 2 ,), since 

'(4 + p 2 )m 2 



g~ 2 (0) = b\n 



(5.8) 



A 2 

Obviously, the maximum value obtained for g 2 p (0), is the one that minimizes (4 + p 2 )ml 
and, at same time, keeps it bigger than A 2 , in order to avoid the appearance of a shifted 
version of the Landau pole. 

Note that the term proportional to p 2 in (|5.4jl is necessary for optimizing the fit for the 
range of momenta q 2 G [0.01 GeV 2 ,2 GeV 2 ], i.e. the region where d(q 2 ) falls down rapidly. 
If we were to set p 2 = in Eq. (j5.4|) we would recover automatically the solution proposed 



m however, such a choice would not correspond to the best possible fit: our numerical 
solution requires bigger values for the coupling g 2 p(<? 2 ), and therefore, p 2 assumes negative 
values, as can be observed on Fig. El where found that p 2 = —3.208. 

It is important to mention that other functional forms for / (q 2 , m 2 (q 2 )) were also tried; 
although some of them could fit d(q 2 ) well, they have been discarded due to the appearance, 
at some scale, of an undesired "bump" in the behavior of g 2 P (<? 2 ). In others words, we have 
required that ^ P {q 2 ) should be a monotonically decreasing function of q 2 . 

In Fig. EH we plot a series of numerical solution obtained by fixing different values for 
d~ l (0). All these solutions have been subjected to the constraint imposed by Eq. (j5.2|) . and 
their corresponding values for a are reported in the inserted legend. Observe that they all 
behave like constants until practically the same scale of 0.01 GeV 2 (the "freezing" plateau), 
then they converge down to a common value at around 1 GeV 2 , beyond which they start to 
develop the known perturbative behavior. 

The running charges, a(q 2 ), for each solution presented in Fig. El are displayed on Fig. El 
Observe that as the value of <i" 1 (0) decreases, the value of the infrared fixed point of the 
running coupling, a(0), increases. Accordingly, from Fig. El and Fig. E3 we can conclude 
that, if smaller values of a were to be favored by QCD, the freezing of the running coupling 
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FIG. 9: Results for d(q 2 ) for various values of d 1 (0) (all in GeV 2 ), together with the respective 
values for a, obtained from the constraint of Eq. ()5.2j) . 

would occur at higher values. It should also be noted that the values of a(0) found here 
tend to be slightly more elevated compared to those of ^l] (for the same value of m%). 

In addition we analyze the dependence of the ratio mo/A on a, extracted from Eq. (|5.fi|) . 
This dependence is shown in Fig. [TTJ corresponding to the cases presented in Fig. We 
observe that as we increase the value of a, namely the sum of the coefficients of the massless 
pole terms appearing in the three gluon vertex, the ratio mo/A grows exponentially as 

^=A 1 exp(j-)+y , (5.9) 

with A x = 0.775, t x = 0.25 and y = 0.436. 

At this point it is natural to ask by how much these results would change if we were to 
turn on again c 2 in Eq. (j5.3|) . We have examined values of c 2 comparable to those used for 
a, i.e. we exclude the possibility that the typically small values of a are a result of a very 
fine-tuned cancellation between large c\ and c 2 , opposite in sign. It turns out that, for all 
previous cases studied, the effect of including non- vanishing c 2 's in Eq. fl5.3|) is less than 3%, 
a fact that justifies their omission. 
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FIG. 10: We plot the correspondent running charge, a(q 2 ), for the gluon propagator presented 
on Fig. EJJ Clearly, we can see as we decrease the values of d _1 (0) the infrared fixed point, a(Q), 
becomes bigger. 

B. Non-canonical solutions 

As we decrease <i -1 (0), it becomes increasingly difficult to fit the numerical solution 
obtained from Eq. ()5.1|) using the expressions given in ()3.17|) . ()5.4|) . and (j5.5|) . In particular, 
for solutions with <i _1 (0) < 0.01 GeV 2 , the discrepancy is relatively large, especially in the 
intermediate region. We will denote such non-canonical solutions by d^s(q 2 ). A typical 
numerical solution is shown by the curve composed by black squares in Fig. It is 
important to emphasize that, although the form of Eq. (j3.17|) is not suited for fitting the entire 
momentum range, one may use it for describing a partial region, i.e. q G [0.2GeV 2 ,A uv ], 
as showed by the line+circle curve on Fig. where we clearly see a quantitative agreement 
between both curves. This observation is related to the discussion following Eq. (|3.25| ). 
allowing us to use the same regularization procedure as in (|5.2j) . 

To see this, let us denote by dc(z) the canonical solution that provides the best fit to a 
given d^(z); then write d^(z) = [dn(z) — dc(z)} + dc(z), and substitute into Eq. (|5.2|) . to 
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FIG. 11: The ratio mo/A as function of the parameter a. The increase is exponential, given by 
Eq.flEU). 

obtain 

d-\0) = ba 

The first term on the RHS of Eq. (j5.10|) is simply the RHS of Eq. (j5.2|) . with d — > dc', it clearly 
converges, since dc is (by construction) a canonical solution. The second term receives an 
appreciable contribution only in the low momenta region, i.e. from to 0.2 GeV 2 , vanishing 
very rapidly in the UV, due to the perfect agreement found between both curves for the 
range [0.2 GeV 2 ,A uv ]. 

Other non-canonical solutions with different <i _1 (0) are plotted on Fig. EE and their 
respective values of a are reported in the legend. Notice that, as the value of d(0) increases, 
the IR plateau becomes narrower. 

Finally, the dependence of <i -1 (0) on a for all cases presented here is shown on Fig. El 
in a log- log plot. The results clearly show a linear behavior for smaller values of <i _1 (0), 
whereas for values of ci^ 1 (0) > 0.01 GeV 2 the growth is exponential. 







d —> dc 


-ba 1 




Jo 



dz z 



d N (z) - d c {z) 
g 2 (z) 



(5.10) 
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FIG. 12: The black square points represent the numerical solution for the RGI quantity d(q 2 ) = 
g 2 A(q 2 ), when we fixed d -1 (0) = 0.001 GeV 2 and a = 7.65 x 10 -3 . The red line + circle curve 
is the fit given by Ea. ()3.17|) . setting p2 = — vr and m 2 , = 0.304 GeV 2 . The small plot is the same 
graph on log- log scale. 

VI. DISCUSSION AND CONCLUSIONS 

In this article we have taken a closer look at various issues relevant to the study of gluon 
mass generation through SD equations. The emphasis of our analysis has focused on the fol- 
lowing points: (i) The gauge-invariant truncation scheme based on the PT-BFM formalism, 
and the possibility it offers in implementing a self-consistent first approximation, omitting 
ghost contributions without compromising the transversality of the gluon self-energy, (ii) 
The necessity of introducing massless poles in the form of the vertex used in the SD equation, 
and the role of the seagull graph in enforcing transversality (iii) A method for regulating 
the resulting seagull-type integrals based on dimensional regularization, and the asymptotic 
properties that the solutions must display in the UV in order for this regularization to 
work, (iv) We have derived a linearized version of the SD equation that originates from the 
one-loop dressed approximation, in the absence of ghost loops, (v) We have introduced a 
phenomenological form for the three-gluon vertex, containing massless poles associated with 
two different tensorial structures, (vi) The resulting equation has been solved numerically 
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FIG. 13: Results for d(q 2 ) fixing different values for d _1 (0) (all in GeV 2 ) in a log- log plot. All 
solutions presented satisfy the condition given by Ea. (|5.1U|) . Their respective values for a are 
described on its legend. 

and two types of qualitatively distinct solutions have been found. 

It should be obvious from the present work that the role of the three-gluon vertex is 

absolutely central; in particular, the question of whether or not it contains massless poles is 

a determining factor when looking for finite solutions. As we have mentioned already, the 

form of the vertex employed here attempts to capture some of the main features, such as 

the impact of the poles, and the possible interplay between poles originating from different 

Lorentz structures (see discussion after (|4.14jl ). Needless to say, an in-depth study of the 

form of the vertex is indispensable for further substantiating the appearance of IR-finite 

solutions. At the moment we only have indications from the study of SD equation and 

lattice simulations (albeit for the three-gluon vertex in the conventional Landau gauge) that 

the IR behavior is indeed singular. In addition, one must explore the possibility of improving 

the gauge-technique inspired Ansatze employed for the vertex, in the spirit of |fi2j], in an 

attempt to correctly incorporate the required asymptotic behavior into the SD equation (see 

I I 

Sec. IIII|) . Given the rich tensorial structure of the three-gluon vertex [68], such a task is 

expected to be technically rather demanding. 
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FIG. 14: d -1 (0) as function of a. 

It would be interesting to further scrutinize the viability of the new class of solutions 
encountered in Sec. El Such solutions are particularly interesting, because they could in 
principle overcome a well-known difficulty related to the breaking of chiral symmetry. Specif- 
ically, the class of solutions displaying freezing, when inserted into standard forms of the 
gap-equation for the quarks, are not able to trigger non-trivial solution, because they do no 
reach high enough values in the IR to overcome the critical coupling [69|. Instead, the new 
type of solutions rises sharply in the deep IR, reaching values that could in principle break 
chiral symmetry, while for intermediate IR momenta it coincides with a canonical massive 



type of solution, that seems to be favored by phenomenology [31| . Of course, it could well 
happen that these solutions are particular to the linearized equation, and do not survive a 
non-linear analysis. 

In recent years the picture that has emerged through the study of SD equations in the 
conventional Landau gauge is characterized by the so-called "ghost dominance" 
In particular, the gluon self-energy has the form A^ u (q 2 ) = A(q 2 )P flu (q), with A(q 2 ) = 
Z(q 2 )/q 2 , and the ghost propagator is D(q 2 ) = —C(q 2 )/q 2 . Assuming that the dressed 
ghost-gluon vertex is finite in the IR, the SD equations yields for q 2 — > 0, C(q 2 ) = Aq~ 2K , 
Z(q 2 ) = Bq AK , where the constants A and B depend on k. The value of k depends on the 
details of the dressing of the gluon-ghost vertex at small momenta; SD and lattice studies 
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seem to restrict it within the range 0.5 < k < 0.6. For the special value k = 0.5 one 
obtains a finite gluon propagator, A(0) = B, whereas the ghost propagator diverges as 
D(q 2 ) — > — A(q 2 )~ ?J l 2 . In addition, by virtue of the identity Z\ = Z g Z\j 2 Z 3 = 1, valid in 
the Landau gauge, where Z 1; Z g , Z 3 and Z 3 are the gluon-ghost vertex, coupling, gluon and 
ghost renormalization constants respectively, one concludes that the product g 2 Z(q 2 )C 2 (q 2 ) 
is RG-invariant, and can be adopted as a definition of the non-perturbative running coupling, 
i.e. a(q 2 ) = (g 2 / An)Z(q 2 )C 2 (q 2 ). Then it is clear that the a(q 2 ) so defined has an IR fixed 
point, a(0) = (g 2 /4:7!-)A 2 B , regardless of the value of k. The actual value of a(0) depends 
on k, through the implicit dependence of A and B on it; for values of 0.5 < k < 0.6, one 
obtains 2.5 < a(0) < 3. Evidently, the (dimensionfull) RG-invariant quantity d(q 2 ) studied 
in this article displays a qualitatively similar behavior to that found for the (dimensionless) 
running coupling in the "ghost-dominance" picture, namely IR finiteness. Therefore, despite 
the difference in the intermediate steps and the terminology employed, the physics captured 
by both pictures appears to be compatible. Notice also that, from the practical point of view, 
the method presented here has the advantage of setting up a SD equation directly for the 
RG-invariant object of interest; thus, the result obtained does not depend on the exactness 
with which a subtle cancellation in the ratio of two quantities, one tending to zero and one to 
infinity, is realized. This last comment may be particularly relevant in the context of lattice 
simulations, where the above ratio is studied numerically; clearly, a small deviation from the 
exact cancellation (due to residual volume dependences, for instance) may lead to serious 
qualitative discrepancies. It is also important to mention that, whereas the confinement 
mechanism within the "ghost-dominance" description is attributed to the divergence of the 



ghost propagator 



7l| . in the picture where the gluon is massive the origin of confinement 



is the condensation of vortices J, . Actually, recent investigations advocate interesting 

I I 

connections between the center- vortex picture and the Gribov-horizon scenario [73j. 

As mentioned in the Introduction, in the PT-BFM scheme the omission of the ghost loop 
does not interfere with gauge-invariance, but it might alter the actual form of the gluon 
self-energy. Therefore, the study of the gluon-ghost system should be eventually considered. 
Of particular relevance in such a study is the nature of the gluon-ghost vertices involved; in 
fact, in the PT-BFM scheme there will be two such vertices: JT U , appearing in Fig. and 
r„ wh,ch will appear in the SD equation for the ghost propagate, The IR behavior of I> 
(in the Landau gauge) is currently under investigation on the lattice |Z4|; it would clearly 
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be important to settle the issue of whether it is divergent or finite. 

Last but not least, the theoretical situation concerning the SD equations within the PT- 
BFM scheme merits further intense scrutiny. First of all, the correspondence between the 
PT and BFM has been established perturbatively to all orders, but no analogous proof exists 
non-perturbatively, i.e. at the level of the SD equations themselves. In this article we have 
assumed that the PT-BFM correspondence persists non-perturbatively. A preliminary study 
in a simplified context (scalar QED) fully corroborates this assumption [75[ ; however, the 
actual realization is highly non-trivial, and its generalization to QCD deserves a thorough 
analysis. In the same context, the consequences of the second crucial ingredient, namely 
the substitution of quantum quantities by background ones inside the loops, are virtually 
unexplored. One must study in detail, preferably in the context of the toy model mentioned 
above, the viability and self-consistency of this procedure. We hope to be able to undertake 
some of these tasks in the near future. 
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APPENDIX: FEYNMAN RULES IN THE BFM 



In this Appendix we list for completeness the Feynman rules appearing in 



{A, a 



r.b 



k k v 

9y.v - (i - io)—^r 



k 2 + ie 



(A-l) 




fj., b 



id 



ab 



k 2 + ie 



gf 



abc 



(Pi - Q) v 9^a + (q~ P2)„ 9ua + (P2 ~ Pl) a 9^ 



(A-2) 



(A-3) 
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a, a 



;i, 6 v, c 



gf 



abe 



pi-q + 7-P2) 9 m +(q-p2- 



SW + (P2 -Pl) a 9tMu 



(A-4) 




-*0 



f aX f X {9iiot9vP — g^/39va) 



+ f adx f xbc (g»»g a p - g^gup) 



(A-5) 




-*0 



1 



/ a x / xc ( gfj.agu/3 — g^pgua + ^-g^g a p 



+ f adx f xbc I g^ v g a p — g m g v f3 — —g^pg, 



[g fn^g &p - g,pg, 



(A-6) 



»iAAAAAA,ft c 
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(A-8) 
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ig 2 f acx f xdh g^ 



(A-9) 
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